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' , We study the renormalization group (RG) running of the neutrino masses and the leptonic mixing 

(— I ' parameters in two different extra-dimensional models, namely, the Universal Extra Dimensions 

(UED) model and a model, where the Standard Model (SM) bosons probe an extra dimension 
^ ! and the SM fermions are confined to a four-dimensional brane. In particular, we derive the beta 

I function for the neutrino mass operator in the UED model. We also rederive the beta function for 

■ the charged-lepton Yukawa coupling, and confirm some of the existing results in the literature. The 
generic features of the RG running of the neutrino parameters within the two models are analyzed 
and, in particular, we observe a power-law behavior for the running. We note that the running of 

' the leptonic mixing angle 9i2 can be sizable, while the running of ^23 and ^13 is always negligible. 

■ In addition, we show that the tri-bimaximal and the bimaximal mixing patterns at a high-energy 
scale are compatible with low-energy experimental data, while a tri-small mixing pattern is not. 
Finally, we perform a numerical scan over the low-energy parameter space to infer the high-energy 
distribution of the parameters. Using this scan, we also demonstrate how the high-energy 9i2 is 
correlated with the smallest neutrino mass and the Majorana phases. 
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I. INTRODUCTION 



With the recent start-up of the Large Hadron Colhder (LHC), experimental physics has 
started the search for the domain beyond the current Standard Model (SM) of particle 
physics. In addition to searching for the Higgs boson, the LHC will supply us with important 
information about the nature of physics above the TeV energy scale. Among the most 
popular high-energy extensions of the SM are extra-dimensional models. The idea that 



spacetime could have more than four dimensions was first proposed by Theodore Kaluza [l] 
and Oskar Klein ^ at the beginning of the twentieth century. In the 1980's, this idea gained 
popularity through the emergence of string theory, and at the end of the 1990's, several 
extra-dimensional models, which could potentially be detected at the next generation of 
high-energy experiments, were proposed 3|-|8|. These models are mainly motivated by the 
fact that they could provide solutions to different problems in the SM, such as the hierarchy 
problem and the lack of a good particle dark matter candidate. 

An interesting aspect of extra dimensions is their impact on the RG running of physical 
parameters. It has been shown that the Kaluza-Klein (KK) towers give rise to an effective 
power-law running of the parameters for energies above the first KK level 91]. Hence, extra 
dimensions may increase the RG running dramatically, resulting in large effects at relatively 
low energy scales. 

One of the rare examples of experimental evidence for physics beyond the SM comes from 
neutrino physics. The observation of neutrino oscillations strongly indicates that neutrinos 
are massive and that lepton flavors are mixed. Since it is not possible to describe neutrino 
masses within a renormalizable framework using the SM particle content only, neutrino 
oscillations imply new physics beyond the SM. The fact that the neutrino masses are bounded 
to be unnaturally small in comparison to the other SM fermion masses, together with the 
possibility that neutrinos could be their own antiparticles, has given rise to extensive model 



building within the neutrino sector. For example, sma^ 
be generated through the so-called seesaw mechanisms 



1 neutrino masses could naturally 



10 



-Il3l]. where neutrinos couple to 



new degrees of freedom. Small neutrino masses could also be generated in certain extra- 



14|-|l7l| , which have the advantage of potentially being observable at the 



dimensional models 
LHC. 

Since the neutrino parameters are measured in low-energy scale experiments, the RG 



running effects should be taken into account properly in studying neutrino mass models at 



high-energy scales 
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- |25| . Therefore, in this paper, we study the RG running of the neutrino 
parameters in extra-dimensional models. We employ an effective description for neutrino 
masses in terms of the dimension-five Weinberg operator, and investigate the RG running of 
the neutrino parameters in two different extra-dimensional models. In particular, we derive 
the beta function for the Weinberg operator in the Universal Extra Dimensions (UED) model 
8|, which has previously not been performed in the literature. In order to determine which 
high-energy values that are consistent with current experimental data within the models, 
we use a Markov Chain Monte Carlo to infer the favored high-energy parameter values from 
the low-energy parameter bounds using our analytical results. 

The rest of the work is organized as follows: In Sec. [Tll we describe the models that we 
have studied, and present the beta functions for the neutrino mass operator in both models. 
In Sec. mil we give the RGEs for the neutrino masses and the leptonic mixing parameters, 
which are obtained from the beta function for the neutrino mass operator. Then, in Sec. IIV[ 
we show our numerical results for the RG running. Finally, in Sec. |Vl we summarize our 
work and state our conclusions. In addition, in Appendix \^ we give the full set of RGEs 
that we have used. 



II. RGES IN EXTRA-DIMENSIONAL THEORIES 



In order to illustrate the features of RGEs in extra-dimensional theories, we consider two 
representative models. First, we investigate the UED model, in which all the SM fields probe 
the extra spatial dimensions. Then, we study a model, in which the SM bosons propagate 
in the bulk, while the SM fermions are localized to a four-dimensional brane jsl. 

In general, the KK towers in an extra-dimensional theory are expected to be cut off at 
some energy scale A in order to keep the theory renormalizable, with the nature of the cutoff 
scale depending on the specific ultraviolet (UV) completion of the model. For the purpose 
of illustration, we take A = 50 TeV for both models in this work. 
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A. Model I: Universal Extra Dimensions 



The UED model is constructed by promoting all the SM fields to a higher-dimensional 
flat spacetime, and hence, all the SM particles acquire towers of KK modes. Here, we 
consider the simplest case with only a single extra spatial dimension, which is assumed to 
be compactifled on an orbifold with radius R. In this framework, KK parity, which is 

deflned as (—1)" for the nth KK level, is conserved after compactiflcation. The mass scale 
of the flrst excited KK level, given by -R~^, is bounded to be larger than approximately 300 



GeV 26|. In this work, we use the value R^^ = 1 TeV. 

In addition to the operators that are renormalizable at the level of the SM, we introduce 
the dimension-flve Weinberg operator responsible for neutrino masses 

- = ^ k {<j)^L') + h.c. , (1) 

where L denotes the lepton doublet flelds, denotes the Higgs doublet, and k is a matrix in 
flavor space. After electroweak symmetry breaking, the neutrino mass matrix is obtained as 

mj, = Kv'^ , (2) 

where u ~ 174 GeV is the vacuum expectation value of the Higgs fleld and n = k/nR. The 
neutrino mass operator deflned in Eq. Q ca n be realized through certain extra-dimensional 



seesaw mechanisms at the cutoff scale A 16|], as well as through the standard seesaw mecha- 
nisms. In the language of an effective theory, Eq. ([1]) is essentially the same for the different 
seesaw models. 

In a flve- dimensional spacetime, there are no chiral fermions. In order to reproduce the 
phenomenology of the SM for energies below R~^, a flve- dimensional Dirac fermion has to be 
introduced for each chiral fermion in the SM, and the KK expansions of these Dirac fermions 
are chosen in such a way that chiral fermions are obtained at the zero-mode level. Hence, 
the number of degrees of freedom in the fermion sector is doubled at the excited KK levels 
in comparison to the SM. In addition, each SM gauge fleld has a flfth component, which 
appears as a real scalar from the four-dimensional point of view. Again, the zero-modes 
of these scalars can be removed by suitable choices of KK expansions, but they appear at 
the excited KK levels. In particular, compared to the SM-like couplings, additional vertices 
involving SM fermions and the flfth components of gauge flelds appear. 



In the UED model, the beta function for the neutrino mass operator can be written as 

16vr2^ = /3SM + ^UED (3) 



where Pf^ denotes the SM beta function [isl, [l9|, l21| 



= -1<Y,%) - - {3g', -2T-X) k, (4) 

where 

T = tr (sY^Yu + 3Y}Yd + Y}Y^ , (5) 

with Yf (for / = n, d, t) denoting the Yukawa couphng matrices of the up-type quarks, 
down-type quarks, and charged leptons, respectively. Here, gi are the gauge couplings and A 
denotes the Higgs self-coupling constant. The second term on the right-hand side of Eq. (|3]), 
comes from the contributions of the excited KK modes. As mentioned above, the KK 
spectrum at the excited levels differs from the SM, and this is reflected in the contributions 
of the KK modes to the beta function. We have calculated with the following result 



UED 



(6) 



where s = [/i//ioJ counts the number of KK levels contributing to the beta function for a 
given energy /i. Here, fio = -R"^ = 1 TeV. For large fi/fio, s is well-approximated by the 
continuous expression fi/fiQ. In evaluating physical parameters from lower to higher energy 
scales, new KK excitations enter the theory at each KK threshold, giving additional quantum 
corrections, and hence, the coefficients of the beta functions are modified depending on how 
many KK modes that are excited. The result is that the RG running of k follows a power-law 
behavior, controlled by s. This is a distinctive feature of RG running in extra-dimensional 
theories. Compared to the RG running behavior in the SM, where the coefficients of the 
beta functions are nearly constant, the power-law behavior results in a significant boost in 
the RG running, which could possibly be tested at near-future experiments. 

The differences in the coefficients for the gauge couplings between Eqs. and (E]) are 
due to the additional Feynman diagrams involving the fifth components of the electroweak 
gauge bosons shown in Fig. [1] The additional factor of 2 in the coefficient for T is due to 
the fact that the chiral fermions are replaced by Dirac fermions at each excited KK level. 




The beta functions for the Yukawa couphngs Yj as well as the gauge couplings are listed in 
Appendix |X1 Note that the beta functions for the Yukawa couplings differ between Refs. [27 1 



and |28j. We have confirmed the computations in Ref. 



UED 



-YYf + 2T 



33 



-9i 



namely, 
15 2' 



(7) 



B. Model II: Fermions on the brane 



In this model, the SM fermions are confined to a four-dimensional brane, while all the 
SM bosons p..obe the bulk g Aga,„. we assume the ext.. di.e„s,o„ to be co.pact.fied on 

an S^/Zo orbifold with radius R. Present collider bounds allow the masses of the lowest KK 

n 

excitations to be as low as 500 GeV [9[. As for model I, we take /io = -R ^ = 1 TeV in this 
work. Thus, the heavy KK modes of the SM particles could be accessible at forthcoming 
collider experiments, and the unification of gauge couplings could also be achieved at a 



low-energy scale 



29 



30|. 



Again, we write the beta function for 



167r' 



dln/i 



(8) 



where the SM beta function is given in Eq. (jl]). The contribution from the excited KK 



modes reads 



91,1311, 



32| 



(9) 



2 V . - 2 

In comparison to the SM beta function, the term proportional to T is missing here. This 
term comes from the fermion loop contributions to the self-energy of the Higgs boson, and 
since the fermions have no KK excitations in this model, there is no such contribution for 
n > 0. Also, the remaining terms in Eq. iQ are larger by a factor of 2, which is due 
to a rescaling of all interactions involvin g K K excitations by a factor v^, coming from the 



canonical normalization of all KK modes 



33| . Note that this factor does not exist in model I. 



III. RUNNING NEUTRINO PARAMETERS 



We proceed our discussion to the RG running of the neutrino parameters. In general, 
one can choose to work in a basis where the charged lepton Yukawa couplin g rn atrix is 
diagonal, i.e., Yi = diag(?/e; 2/^? ^r)- In this basis, the leptonic mixing matrix 3J, |35| stems 
from the diagonalization of the neutrino mass matrix, i.e.. 



U'^m,yU* = Dy = diag(mi, m2, ms) 



(10) 



with rrii being the neutrino masses. Inserting Eq. f lTOj) into the beta function for k, one 
obtains the evolution for the leptonic mixing matrix due to the KK modes as 

dU 



dt 



U = 2sUX , 



with t = ln/i/167r^, and 



ImX,, = -^r/Q^Im (u^^YeU 



(11) 

(12) 
(13) 



2 - \ / ij 

where 77 = 1 in model I, and 77 = 2 in model II, and the indices i and j run over 1, 2, and 



3. The factors Qj are defined as C, 



[mi 



mi)/{mj — rrii). Note that the charged lepton 



mass spectrum is strongly hierarchical, i.e., rrie ^ ^ m,-, which allows us to make a 
reasonable approximation by ignoring the electron and muon Yukawa couplings in Eqs. f lT2|) 
and (ITS]) . In what follows, we will assume = diag(0, 0, t/r) for simplicity. Furthermore, the 
C factors play a key role in the RG running of the neutrino mixing angles, since in the case 
of a nearly degenerate neutrino mass spectrum, i.e., mi ^ m2 — m^, Qj ^ 1, and therefore, 
the RG running effects will be enhanced dramatically. If the neutrino mass spectrum is 
hierarchical, e.g., mi <^ m2 <^ ma, Qj ~ 1 holds and the RG running effects on the leptonic 
mixing matrix are not observable. In the following analytical analysis, we will assume a 
nearly degenerate neutrino mass spectrum. 

In order to figure out the RG running behaviors of the leptonic mixing parameters, we 
employ the standard parametrization, in which U is parametrized by three mixing angles 
and three CP-violating phases as 



U 



C12C13 



S12C13 



-S12C23 - £125235136"^ C12C23 - 5123235136"^ 



Sl3e 
S23C13 



\ 



S12S23 - Ci2C23Sl3e' 



-C12S23 - Si2C23Sl3e'" C23C13 



v 



V 



(14) 



with Cij = cos dij and Sij = sin 6ij {ij = 12, 13, 23). Combining Eqs. ( ITTj) and ( !T4|) . we arrive 
at the running for the leptonic mixing angles 



e 



12 



VsCl2Sl2Cl2sl3cl_^yr , 



^23 - T7^^Cl3S23C23 (^igcj + cl^cl) y. 



'23 
^13 



2 



2 

2'7'5Cl3Sl2Cl2S23C23 (c^Q+a - CpQ+p) y. 



2 



(15) 
(16) 
(17) 



where = cosx. Here, we have taken (13 ^ (23 and neglected the smallest mixing angle 
6*13 to a reasonably good approximation. In the limit of four dimensions, i.e., s = t] = 1, 



the RGEs for 9ij in the SM are reproduced [36|, l37|. One observes that, compared to the 



RG running of the mixing angles in four-dimensional theories, an additional enhancement 
factor s enters the RGEs, which may lead to significant RG running effects. However, as 
the beta functions for the mixing angles are proportional to the four- dimensional ones, some 
qualitative features (e.g., the dependences of the running mixing angles on the neutrino 
mass hierarchy and the Majorana CP-violating phases) remain the same, and were already 



known in four- dimensional theories 



18 



-|25|. 



For the RG running of the leptonic mixing angles, one of the key features is that 612 
increases with increasing energy scale, independently of the neutrino mass hierarchy and 
the CP-violating phases. Consequently, a small 612 at the high-energy scale is disfavored. 



On the other hand, the tri-bimaximal (si2 = l/"\/3, S23 = l/v^? -^is = 38l-|40| ) and the 
bimaximal (si2 = S23 = l/V^, S13 = 0) mixing patterns could be natural candidates for 
flavor symmetries at some high-energy scale. As for 623 and ^13, the RG running effects are 
milder, since their RG running is boosted by C13, which is much smaller than (12- 

For the sake of completeness, we also give the analytical RGEs for the neutrino masses 

"^1 - -^Vsmisl^sl^yl + rrixa^ , (18) 

rn2 ^ ~'nsm2sl^cl2yl + m2a,, , (19) 
3 

rn-s ^ --rjsm-icl^yl + m-iU^ , (20) 

where is flavor universal, and can be found in Appendix [A] Similarly to the RG running of 
the mixing angles, the power-law factor s enhances the RG corrections. However, compared 
to the RGEs for the mixing angles, there is no enhancement factor C,. Furthermore, the flavor 
non-trivial parts in Eqs. (lT8l)-(l20l) are suppressed by the charged- lepton Yukawa coupling yr 
in comparison to In particular, if one ignores the terms proportional to y^, and obtain a 
very simple form for Eqs. (fT8l) - (!20l) as rrij ~ Tnia,^. The solution to this equation is roughly 
estimated by 

mi(A) / A 



. ^ - , , • (21) 

mi{Mz) \Mz) ^ ' 

Therefore, the RG running of the neutrino masses is only sensitive to a^? independently of 
the neutrino mass spectrum and the mixing parameters. 

In the following numerical analysis, we will mainly concentrate on the RG corrections to 
the flavor structure of the leptonic mixing matrix in the two above models. 

IV. NUMERICAL ANALYSIS 

In our numerical computations, we make use of the full sets of RGEs without any ap- 
proximations. The input values for the neutrino parameters and SM observables at the 



yU. = Mz scale are taken from Refs. 4l|, |42]. In addition, for the Higgs mass, we use the 



representative value ntH = 140 GeV. Direct information on the absolute neutrino mass scale 



60 



45 



.30 



15 



10 



60 



45 



.30 



15 



10 



Normal Hierarchy 



Normal Hierarchy 



10 10 

/i [GeV] 



12 

^3 



10" 10' 

fi [GeV] 



60 



45 



.30 



15 



10' 



60 



45 



.30 



15 



10 



Inverted Hierarchy 



Inverted Hierarchy 



10 10 

/i [GeV] 



10" 10' 

^l [GeV] 



FIG. 2. The RG evolution of the three leptonic mixing angles as functions of the energy scale 
from Mz to A for the normal mass hierarchy (left column) and the inverted mass hierarchy (right 
column) in model I (upper row) and model II (lower row). The solid, dashed, and dotted curves 
correspond to the mass of the lightest neutrino to be 0, 0.2 eV, and 0.5 eV, respectively. Since 
there is no strong experimental evidence for a non- vanishing ^13, we take ^13 = as the input value 
at the A4z scale. In addition, all the CP-violating phases are taken to be zero. 



can be derived from tritium beta decay experiments 43|, |4J], i.e., < 2.3 eV (at 95 % 
C.L.). Indirect constraints from the CMB data of the WMAP experiment and the large 
scale structure surveys also lead to an upper limit on the sum of neutrino masses. In our 
numerical studies, we conservatively take nii < 0.5 eV. 

In Fig. [21 we show the RG evolution of the neutrino mixing angles. In good agreement 
with our analytical results, 612 is the mixing angle, which is the most sensitive to the RG 
corrections, and it increases with increasing energy in both the normal hierarchy (mi < 
7712 < '^3) and inverted hierarchy (7713 < mi < m2) cases. In contrast, 623 and ^13 are rather 
stable under the RG running, which reflects the fact that there is no strong enhancement 
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FIG. 3. The dependence of 612 on mi at A = 50 TeV in model I. The band between the two curves 
corresponds to the allowed parameter space at 90 % posterior probability. Here, a normal mass 
hierarchy is assumed. 



factor C12 in their RGEs. Furthermore, the running direction of 623 depends on the mass 
hierarchy, namely, an increasing (decreasing) 6*23 is obtained in the case of normal mass 
hierarchy (inverted mass hierarchy). This feature can be understood from the sign of C13, 
which is positive in the normal hierarchy case, and negative in the inverted hierarchy case. 

Note that the RG running effects on the mixing angles in model I are generally larger 
than those in model II, although the coefficients of the beta functions in model II are twice 
as large due to 77. This can be understood from the different RG running behavior of t/r in 
the two models. Concretely, in model I, the flavor independent parts of Py^ contain the trace 
of charged- fermion Yukawa couplings (c.f., the T term in Eq. (lAlSP ). whereas such a contri- 
bution does not exist in model II, due to the absence of fermion KK excitations. Therefore, 
t/r receives larger RG corrections in the UED model and a more sizable yr could naturally 
be expected at higher scales, which eventually leads to more signiflcant RG corrections to 
6ij in the UED model. 

In Fig. [31 we present the dependence of the high-energy 612 on mi in the normal mass 



hierarchy case. In the numerical calculations, we use the MonteCUBES software 45|] as a 
basis for a Markov Chain Monte Carlo (MCMC), and generate 5 ■ 10^ samples per value of 
the mass mi, while the neutrino oscillation parameters have priors corresponding to present 
bounds. The band between the two curves corresponds to the allowed parameter space at 
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FIG. 4. The RG evolution of the ratio mi{K)/mi{Mz) on mi for model I (solid curve) and model 
II (dashed curve). Here, we use mi = 0.2 eV and a normal mass hierarchy. 

90 % posterior probabilitjll]. One can observe that remarkable changes in 6*12 can be achieved 
at the higher energy scales if mi is sufficiently large. For a proper choice of mi, the tri- 
bimaximal or the bi-maximal mixing pattern can be easily achieved at higher-energy scales. 
However, in contrast to the conclusion obtained in Ref. 32j, large leptonic mixing angles at 



low-energy scales cannot be generated from very small mixing angles near the cut-off scale. 

This is due to the naive but unrealistic two-flavor picture employed in Ref. 32[. Explicitly, 

in the two-flavor framework, a small mixing angle can be obtained at the cut-off scale for 

m^/Am^ ~ 10'^. Thus, when we are working in the standard three-flavor framework, the 

RG corrections to 6*13 and 62^ are related to {m^ + m2)^/Am32, which can be at most at 

order 10^ for m^ < 0.5 eV, and hence not sufficiently large to result in visible effects. As for 

^12, the RG running is related to the ratio (m2 + mi)^/Am2i, which could be large enough 

to give significant corrections for a nearly degenerate neutrino mass spectrum. However, as 

we have mentioned, the RG running always leads to a larger value of 6'i2 at the cut-off scale, 

independently of the choices of neutrino mass hierarchy. Therefore, the tri-small mixing 

pattern is not compatible with the models under consideration. We further remark that this 

feature is generic and independent of the choices of other physical parameters. 

For the sake of completeness, we illustrate the running of the neutrino masses in Fig. HI As 

expected, the RG running behavior is universal for the three neutrino masses within a specific 

^ The 90 % posterior probability region is the smallest region containing 90 % of the posterior probability 

distribution, i.e., given the priors, there is a 90 % probability that the parameters are within this region. 

It is to Bayesian statistics what confidence regions is to frequentist statistics. 
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model. For model I, the RG corrections could increase the neutrino masses dramatically, 
whereas in model II, the RG corrections are not significant for the neutrino masses. This is 
consistent with our analytical result that is larger in model I, due to the contributions 
to T from the top quark Yukawa coupling. We have also checked that this conclusion does 
not change for different hierarchies of the neutrino masses. 

Finally, since the RG running of 612 is very sensitive to the Majorana CP- violating phases, 
it is of interest to observe the correlations between the CP-violating phases and the RG 
running of the leptonic mixing angles. To this end, we show in Fig. [5] the correlations 
between 6*12 at the cutoff scale and the phase difference p — a, with 5 ■ 10^ samples per value 
of p — cr. The bands in the plot indicate the favored parameter spaces (at 90 % posterior 




0.5 1 1.5 2 0.5 1 1.5 2 



lp-al/7i lp-Ol/7I 

FIG. 5. The 90 % posterior probability regions for 012(A) as depending on the Majorana CP- 
violating phase difference p — a in model I (left panel) and model II (right panel). The bands 
between red, green, and blue curves correspond to the allowed parameter space for mi = 0, 0.2 eV, 
and 0.5 eV, respectively. 

probability) for mi = (red), 0.2 eV (green), and 0.5 eV (blue), respectively. As expected, 
no visible RG running effects are observed for a hierarchical neutrino spectrum. As for the 
degenerate neutrino spectrum, sizable RG corrections may exist depending on the difference 
between two Majorana CP-violating phases. A peak in 6'i2(A) appears around p — cr = nvr, 
while the RG corrections are damped if p — cr = (n -|- l)7r/2, with n being an integer. 
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V. SUMMARY AND CONCLUSION 



In this work, we have studied the RG running of neutrino parameters in two representative 
extra-dimensional models extended by an effective neutrino mass operator. In particular, we 
have derived the full set of RGEs for the Yukawa coupling matrices and the neutrino mass 
operator. Both analytical and numerical analyzes of RG corrections to the neutrino mixing 
angles and masses have been performed according to our RGEs. We have found that, due to 
the power-law behavior and a sizable enhancement factor, 612 is the most sensitive mixing 
angle to the RG running, especially in the nearly degenerate limit, e.g., 623 = vr/4 can be 
accommodated at high energy scales if mi ~ 0.5 eV. In contrast, the RG corrections to 
623 and 613 are negligible in both of the models. In addition, the RG running effects on 6*12 
in model I are generally larger than those in model II, due to the charged-lepton Yukawa 
couplings. Most interestingly, 6*12 does not decrease with increasing energy scale, regardless 
of the choice of the neutrino mass hierarchy and the CP-violating phases. Therefore, mixing 
patterns with small 6*12 at the cutoff scale are not compatible with low-energy experiments, 
whereas the tri-bimaximal and bi-maximal patterns turn out to be favorable, depending on 
the specific choice of model parameters. We have also presented the connection between the 
running 6*12 and the Majorana CP-violating phases, which indicates that the RG correction 
is damped if the difference between the Majorana phases is close to (2n-|-l)7r/2. We conclude 
that quantum corrections should not be neglected in studies of extra-dimensional neutrino 
mass models. Our results allow to carry out an integrated investigation of fermion masses 
and flavor mixing in the framework of extra dimensions. 

In the current work, we have only considered theories with one spatial extra dimension. 
In the cases of models with more than one extra dimension, the RG corrections are generally 
more substantial. The conclusions of our study depend on the compactification scheme, and 
other choices of boundary conditions would lead to different phenomena, which are however 
beyond the scope of the current work. 
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Appendix A: Full set of one-loop RGEs for the extra-dimensional models 



The one-loop RGEs for the Yukawa couphng matrices Yf (/ = u, d, i) and the neutrino 
mass operator k, can be expressed in a general form as 



167r' 



dF, 



dln/i 



dln/i 

where the SM beta functions are [18 

'3 



(Al) 
(A2) 



19 



2l| 



SM 



Pi 



17 2 9 2 



gi - 8gi + T 



SM 



(A - 3gl + 2T) «: 



15 2 ^2 



-9i 



3 
-/ 
2 



4 

3 / X \ 7" 



(A3) 
(A4) 
(A5) 
(A6) 



with T = tr (^Y^Yu + 3FjF(i + Y^Y^ . The contributions from KK excitations are given by 



^VsYX - IvsYjY, , 
-^r^sY^Y^ + ^T^sYjY,., 
IvsY^Y, , 
—VsY.^Ye, 



(A7) 
(A8) 
(A9) 
(AlO) 



where rj = 1 for model I and t] = 2 for model II. Here, we have defined the scale parameter 
s = [fi/fio\. The flavor diagonal coefficients a's for model I read 



au = Tjs 
Oid = vs 
an = rjs 



101 2 15 2 28 , ^ 
—g,--g2--g, + 2T 



17 



15 



28 



--g,--g,--g, + 2T 



Ok = Tjs [ -^gf - ^gl + 4T + A 



(All) 
(A12) 
(A13) 
(A14) 



16 



whereas for model II, we have 





= rjs 


( 17 2 
V 12 


9 2 

4 


Old 


= rjs 


( 5 2 
1,-12^^ - 


9 2 


at 


= TjS 




9 2 




= rjs 


(A - ?>gl) . 





^72 - 8^73 



For the sake of completeness, we also present the RGEs for the gauge couplings 



167r 



2 



dln/i 



(A15) 

(A16) 

(A17) 
(A18) 

(A19) 



where {bf^ ,bf^ ,hf^) = (41/6,-19/6,-7), while (61,62,^3) = (27/2,7/6,-5/2) in model I 
and {bi, 62, ^3) = (1/6, —41/6, —21/2) in model II. Furthermore, we also need the running of 
the Higgs self-coupling, 

dA 



167r' 



where the SM contribution reads 46 1 



dln/i 



(A20) 



oSM 

Pa 



- A {3gf + 9gl) + (^^gt + Sg^gl + ^gt 



+4AT - 8tr 

In addition, the extra-dimensional contributions are 



(A21) 



/3a = QvsX^ - VsX {^gl + 9gl) + r/s {2gt + Aglgl + Qgt) 



-SrjsXT — 16?7str 



3 {Y^Y^f + 3 (YjY.y + (y^Y, 



for model I, and 



= Qvs>? - VsX {3gl + 9gl) + r]s ( ^gt + 3glgl + "^-g^ 



(A22) 



(A23) 



for model II. 



[1] T. Kaluza, Sitzungsber. Preuss. Akad. Wiss. Berlin (Math. Phys. ) 1921, 966 (1921). 
[2] O. Klein, Z. Phys. 37, 895 (1926). 

17 



[3] 
[4] 
[5] 

[6] 
[7] 
[8] 

[9] 

[10] 
[11] 

[12] 
[13] 

[14] 
[15] 

[16] 

[17] 

[18] 
[19] 

[20] 

[21] 



I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos, and G. Dvali, |Phys.Lett.| B436, 257 (1998), 
arXiv:hep-ph/9804398 [hep-ph ][ 



N. Arkani-Hamed, S. Dimopoulos, and G. R. Dvali, Phys. Lett. B429, 263 (1998), 



hep-ph/9803315 



N. Arkani-Hamed, S. Dimopoulos, and G. R. Dvali, Phys. Rev. D59, 086004 (1999), 



hep-ph/9807344 



L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999), hep-th/9906064 



L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999), hep-ph/9905221 



T. Appelquist, H.-C. Cheng, and B. A. Dobrescu, |Phys. Rev.| D64, 035002 (2001), 



arXiv:hep-ph/0012100 



K. R. Dienes, E. Dudas, and T. Gherghetta, |Nucl. Phys.| B537, 47 (1999) 



arXiv:hep-ph/9806292 



P. Minkowski, Phys. Lett. B67, 421 (1977). 

T. Yanagida, in Proc. Workshop on the baryon number of the Universe and unified theories, 

edited by O. Sawada and A. Sugamoto (1979) p. 95. 

R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980). 

M. Gell-Mann, P. Ramond, and R. Slansky, in Supergravity, edited by P. van Nieuwenhuizen 
and D. Preedman (1979) p. 315. 



J. M. Prere, G. Moreau, and E. Nezri, [Phys7Rev:| D69, 033003 (2004), |arXiv:hep-ph/0309218 



N. Haba, S. Matsumoto, and K. Yoshioka, |Phys. Lett. I B677, 291 (2009), 



arXiv:0901.4596 [hep-ph 



M. Blennow, H. Melbeus, T. Ohlsson, and H. Zhang, Phys. Rev. D82, 045023 (2010) 



arXiv: 1003.0669 [hep-ph 



T. Saito, M. Asano, K. Fuju, N. Haba, S. Matsumoto, et al, Phys. Rev. D82, 093004 (2010) 



arXiv: 1008.2257 [hep-ph 



P. H. Chankowski and Z. Pluciennik, |Phys. Lett. B316, 312 (1993), arXiv:hep-ph/9306333 



K. S. Babu, C. N. Leung, and J. T. Pantaleone, |Phys. Lett.| B319, 191 (1993) 



arXiv:hep-ph/9309223 



S. Antusch, M. Drees, J. Kersten, M. Lindner, and M. Ratz, Phys. Lett. B519, 238 (2001) 



arXiv:hep-ph/0108005 



S. Antusch, M. Drees, J. Kersten, M. Lindner, and M. Ratz, Phys. Lett. B525, 130 (2002) 



18 



arXiv:hep-ph/0110366 



W. Chao and H. Zhang, iPh^. Rev. ID 75, 033003 (2007), [a^v:hep-ph/06lT323 



M. A. Schmidt, [Phys. Rev.| D76, 073010 (2007), |arXiv:0705.3841 [hep-phj 



J. Chakrabortty, A. Dighe, S. Goswami, and S. Ray, |Nucl. Phys.| B820, 116 (2009) 



arXiv:0812.2776 [hep-ph; 



J. Bergstrom, M. Mahnsky, T. Ohlsson, and H. Zhang, |Phys. Rev.| D81, 116006 (2010) 



larXiv: 1004.4628 [hep-ph 



D. Hooper and S. Profumo, |Phys. Rept., 453, 29 (2007), arXiv:hep-ph/0701197 [hep-ph 



G. Bhattacharyya, A. Datta, S. K. Majee, and A. Raychaudhuri, |Nucl. Phys.| B760, 117 



(2007), arXiv:hep-ph/0608208 



A. S. Cornell and L.-X. Liu(2010), arXiv:1010.5522 [hep-ph 



I. Antoniadis, |Phys. Lett.|B246, 377 (1990). 



K. R. Dienes, E. Dudas, and T. Gherghetta, |Phys. Lett.| B436, 55 (1998) 



arXiv:hep-ph/9803466 [hep-ph 



S. A. Abel and S. F. King, |Phys. Rev:] D59, 095010 (1999), [a^v:hep-ph/9809467 



G. Bhattacharyya, S. Goswami, and A. Raychaudhuri, Phys. Rev.| D66, 033008 (2002) 



arXiv:hep-ph/0202147 



T. G. Rizzo and J. D. Wells, |Pgs[Re^D61, 016007 (2000), [aIxrv:hep-ph/9906234] 



Z. Maki, M. Nakagawa, and S. Sakata, [Prog. Theor. Phys] 28, 870 (1962). 
B. Pontecorvo, Sov. Phys. JETP 26, 984 (1968). 

S. Antusch, J. Kersten, M. Lindner, M. Ratz, and M. A. Schmidt, IJHEPI 03, 024 (2005), 



arXiv:hep-ph/0501272 



J.-w. Mei, |Phys. Rev^ 'DTl, 073012 (2005), [aIOCrv:hep-ph/05020T5 



P. F. Harrison, D. H. Perkins, and W. G. Scott, |Phys. Lett.] B530, 167 (2002) 



arXiv:hep-ph/0202074 



P. F. Harrison and W. G. Scott, |Phys. LetEj B535, 163 (2002), arXiv:hep-ph/0203209 



Z.-z. Xing, |Phys. Lett] B533, 85 (2002), arXiv:hep-ph/0204049 



T. Schwetz, M. A. Tortola, and J. W. F. Valle, [New J. PhysT] 10, 113011 (2008) 



arXiv:0808.2016 [hep-ph 



Z.-z. Xing, H. Zhang, and S. Zhou, Phys. Rev.|D77, 113016 (2008), arXiv:0712.1419 [hep-ph; 



V. Lobashev, Nucl. Phys. A719, 153 (2003). 



19 



[44] K. Eitel, |Nucl. Phys. Proc. Suppl.| 143, 197 (2005). 

[45] M. Blennow and E. Fernandez-Martinez, |Comput. Phys. Commun.| 181, 227 (2010), 



arXiv:0903.3985 [hep-ph] 



[46] T. P. Cheng, E. Eichten, and L.-F. Li, |Phys. Revl D9, 2259 (1974). 



20 



